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1 Introduction 



The Euler equations appear in hydrodynamics as scaling limits for the dynamics of the conserved 
quantities in a fluid. So far the only rigorous example of a microscopic derivation starting from 
Newton's laws is that of the one-dimensional hard rod fluid, [BDS|, |d), pi] , [BS]. In this degen- 



erate situation the number of particles with a given momentum is locally conserved. Therefore the 
hydrodynamic field n v (x, t) is just counting the number of particles at time t with space coordinate 
x and with momenta equal to v (we have put the masses equal to one). When the rod length tends 
to zero we deal with a one-dimensional ideal gas and, after the appropriate (Eulerian) scaling, the 
associated density field converges to the solution of 

d t f{x,v;t)+vd x f(x,v;t) = (1.1) 

(for certain initial conditions). One can then go further and study the fluctuations around this 
macroscopic equation. They are Gaussian with a covariance matrix containing information about 
the transport coefficients which are finite and non-zero for the hard rod fluid. 

In the present paper we investigate what remains of this in the quantum case. Then also here we 
pick out the simplest possible quantum dynamics and we simply ask how the corresponding Euler 
equations and the fluctuations around it can be derived. We see this as a modest step in the rigorous 
study of what could be called quantum hydrodynamics. Clearly, a mathematical derivation ab initio 
of the macroscopic equations of so-called quantum liquids is beyond reach. We may however learn 
something about the conceptual set-up to start such a project from a rigorous study in the most 



simple cases. We refer to [|S1|], [32 for an overview and some solutions to some of the questions 
in the quantum domain of hydrodynamics. In all cases the underlying philosophy is however not 
fundamentally different from the classical situation. Conserved quantities vary on a much larger 
time scale than the others and by a law of large numbers the macroscopic equations appear as closed 
equations governing the motion of the rescaled conserved quantities. 

Here we study quasi-free systems. That means that the microscopic dynamics is in some sense 
linear (as for classical Gaussian systems) and the equations are determined by the one-particle 
motion. Therefore we do not have to deal with the most important problem of time-scale separation 
(the Boltzmann-Gibbs principle) to see how the rescaling effectively truncates the hierarchy of kinetic 



equations. In this situation two papers [S2fl and [S3 1 by Herbert Spohn show how to get the program 



started. In the next Section we summarize some of the ideas. In Sections 3 and 4 we discuss the 



derivation of the Euler equation (similar to (1.1)) for free systems subject to a slowly varying external 
potential. It is interesting to see that this can be applied (in Section 5) in the case of hard core 
bosons on the one-dimensional lattice. Section 6 deals with the non-equilibrium fluctuations. We 
compute the covariance of the rescaled density field. As expected, only at that moment do we start 
seeing a difference between the quantum statistics. Finally, in Section 7, we turn to the case of 
non-interacting quantum particles in a periodic potential. 



2 Some Preliminaries 

The first step in doing anything related to quantum hydrodynamics is probably asking for what 
quantum systems we can get a good idea of what are the conserved quantities and what is the 
structure of the equilibrium states. The answer is disappointing: compared to the situation in so- 
called classical statistical mechanics, there is very little we know about truly interacting quantum 
systems. It is therefore not unreasonable to turn first to quasi-free systems where already on the 
microscopic scale the dynamics is governed by a one-particle Hamiltonian. A quasi-free dynamics is a 
quantum dynamics generated by a Hamiltonian H which is quadratic in the creation and annihilation 
operators a*(x), a(x) of a Bose or Fermi field. Formally, 

H = J dxdy h e (x,y)a*(x)a(y). (2.1) 

As usual with a varying number of particles, this should be thought of as acting on Fock space 
with single particle space L 2 (R d , dx). The quadratic form is specified by a one-particle Hamiltonian 
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h s which in the physically more interesting situations has the form 



-[p + eA(ex)] 2 



U(x) + V(ex) 



(2.2) 



where U is a periodic potential, V is an external potential and A is the vector potential. The 
parameter e > specifies the scale over which the external potentials are varying and will be our 
scaling parameter in what follows. In the next Section we start with the case where U = but 



allowing a bit more general form of h s (x,y) than obtained fr om (2.2). It is not too difficult to see 
that, when dealing with a quantum dynamics generated by (2.1), the evolution of the correlation 
functions can be specified in terms of the single particle evolution generated by h e . Therefore, the 
kinetic equations are closed even on the microscopic level which is an enormous advantage for doing 
hydrodynamics. 



Notice however that the form (2.1) need not be restricted to the microscopic domain (where 



(2.1)) is rather simplistic); there are plenty of quasi- free models in mesoscopic physics where (2.1 
appears as effective action or as describing the evolution of quasi-particles. E.g. the central question 



in Bogoliubov's work on superfluidity was to see how to arrive at something like (2.1) (with some 
very particular properties) starting from a realistic interaction between He 4 -particles, AVZ[. 
This is just a question of equilibrium statistical mechanics. Once arrived at (2.1) with an energy 
spectrum E{k) ~ \k\ for small |fc|, Landau's theory (as e.g. in |LP| ) takes over and completes the 
dynamical picture of superfluidity. 

A quasi-free state is a state on the CAR or CCR algebra (the algebra's generated by the identity 
and the a*(x), a(x) satisfying the Canonical (Anti) Commutation Relations) for which all c orre lation 
functions (and hence the state itself) are determined by the two-point functions, see e.g. [BR]. For 
gauge invariant quasi-free states u> we must only specify the (v(a*(x)a(y)), which, in turn, are given 
in terms of a self-adjoint operator on the one-particle space. 

Because of the dependence on the one-particle evolution in our systems (as explained abo ve) it 
is crucial to understand the scaling of the long time beha vior of a quantum particle subject to ( |2.2| ). 
To keep it simple, let us for a moment put A — in (2.2). Consider the position operator r(t) in the 
Heisenberg picture and define its rescaling as 



r £ (t) = er(e~ 1 t) 



(2.3) 



Rescaling the momentum operator as p e (t) = p(e t) we see that the commutator [r £ ,p e ] vanishes 
as e goes to zero. Let us therefore consider the semiclassical equations of motion, cf. [GMMP|, [|R|], 
[MMP|, (Si. Writing E n (k) (n = 1, 2, . . . and k 6 the Brillouin zone), for the (Bloch-)eigcnvalues of 



the periodic part h° = ^p + U(x) (with periodic boundary conditions), for each band the semiclas- 
sical equations (K = 1) look as follows: 



d t r = d k E n {k), 
d t k = -d r V(r). 



(2.4) 
(2.5) 



If at time t = the particle has a probability distribution / ra (r, k), then, via (2.4-2.5), the distribution 
at time t > is given by / ra (r, k; t) = f n (r(—t), k(—t)) and solves 



d t f n (r,k;t) +d k E n (k)d r f n (r,k;t) = d r V(r)d k f n (r,k;t). 



(2.6) 



We can therefore expect that (|]^) gives the correct Euler equation for the considered quasi-free 
system but we still need to understand what to use for distribution function f n (f, k;t). 

The next question is thus to see what is the analogue of classical quantities like n v (x, t) or f(x, v; t) 
(1.1). Let us therefore briefly recall the notion of Wigner distribution function (cf. fwf[) . 



appearing m 



A recent and to this work very relevant mathematical survey is contained in Section 1 of ]GMMP ] 
If ip( r ) is t ne wave function of a quantum system, then we call (h = 1) 



F(r,k)= - 

TT 



drj e 2lk ' v ip(r + r])tp(r - rj), 



R-' 



(2.7) 
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the distribution function of the simultaneous values of the coordinates r and momenta k. Even 
though F is not positive in general (but it is real valued), this is traditionally justified by the 
following properties: 



I^MI 2 = / dkF(r,k), 
|^(fc)| 2 = / drF{r,k). 



(2.8) 
(2.9) 



By~we denote Fourier transformation. 

We wish to do exactly the same thing for our quasi-free systems. We will come back to the case 
of a periodic potential in Section 7 but for the moment we now also put U — 0. Now we have, 
for every e > small, a state uj e and we must rescale space and time by e . Let a* denote the 
microscopic time evolution to be generated by a quasi-free Hamiltonian and write the time-evolved 
state as lo £ o a s -i t — wf_ lt . Then we define (cf . p2| ) , 



f(r,k;t) 




dr ] e 2ir >- k uj e e - H [a*(E- 1 r + rj)a(£- l r - ij)] 



(2.10) 



as the macroscopic (scale e~ l ) one-particle distribution function of the system. Clearly, clustering 
conditions on the state uj e will be required to allow the convergence of the integral (sum) over rj. 
From now on the assumption stands that this clustering holds allowing a well defined (2.10). Even 
though we write an integral in (2.10) we prefer to avoid here technicalities related to working in the 
continuum. We therefore mostly think of lattice systems (in which case we really have a discrete 
sum, z r should be replaced by its integral part and the wave vector k is in the first Brillouin 
zone). Expressions like a*(x)a(y) should be understood in the distributional sense. Also later we 
will use the int egr al sign as a common symbol even when considering discrete systems. The analogous 



properties to (2.8)-(2.9) apply. As an example, note that if we are interested say in the density and 



0) we have a product state with oj £ [a*{e r + rf)a{e r — rj)] — 8{rj)p(r), then 
is the initial particle density. 



if initially (at t = 
/ e (r,fc;0) = £pW 

The main question is now to investigate the limiti ng b ehavior (as e goes to zero) of ( 2.1C ). That 
is, to derive what corresponds to the Euler equation (1.1). Afterwards we look for the fluctuations 
around this limiting behavior. 

It should be clear by now that we do not identify the problem of deriving 'quantum' Euler equa- 
tions or of studying quantum hydrodynamics with that of 'quantizing' the classical hydrodynamic 
equations. The macroscopic equations remain just 'ordinary' partial differential equations for con- 
served quantities (and not for operators). Still, the quantum nature of the underlying system may in 
principle have a non-trivial effect on these equations (as for quantum liquids, superfluidity etc.); if 
not on the macroscopic equations themselves then on the fluctuations (with corresponding transport 
properties) around it. 

Scaling limits for Wigner functions (measures) are studied in the mathematics literature (see for 
instance the recent work by [ GMMP | ) . The starting point there is the Schrodinger equation 



dt^ie-'x, £-H) + ih e ip e {e- x x, £~H) = (2.11) 

with ip E (x,t = 0) = ^ e (:r)n. h E is the same type of Hamilton operator as in our quasi- free systems. 
In fact, we shall use their result to show the existence of f(r,k;t), the limit of (2.10) as e goes 
to zero. Moreover, to a large extent hydrodynamics of quasi-free quantum systems reduces to the 
homogenization limits considered in |GMMF |. 



1 Regarding notation, we identify the u £ (ex,et) of [GMMP with our ^(x,t) 



3 



3 Euler Equation for Quasi-Free Systems 



We consider a system of Bosons or of Fermions with formal Hamiltonian 

H = I dxdy h e (x,y)a* (x)a(y). 



As usual, the field of creation and annihilation operators is denoted by a*(x), a(x). The kernel h e (x, y) 
corresponds to a one-particle Hamiltonian. We always assume here that h £ (x, y) = h(e(x+y)/2, y—x) 
with an appropriate decay condition, e.g. sup^ J dy\h e (x 1 y)\ < oo for Fermions, to generate an 
infinite volume quasi- free time evolution further denoted by at- Or put differently, the matrix 
elements h s (x — + |) = h(ex,y) may be thought of as hopping rates (from place x — U to 
x + | by a distance y) varying slowly with x. The time t is real. Equivalently, we may write the 

Hamiltonian H via the Fourier transform d(p) = (t^:)^ 2 J due lpw a{u), 

H= (^-\ JdkJdpJ due l(k - p)u E(eu,(k+p)/2)d*(k)d(p), 

where the energy spectrum, E, is defined by: E(r,k) = J dv e~ lkv h(r,v). Clearly, in the limit of 
infinite scale separation, e to zero, the microscopic time evolution is translation invariant with (true) 
energy spectrum E{k) = E(0,k). Yet, as we will see shortly, the (macroscopic) Euler equation 
remembers the dependence of E(r, k) on r. Finally, the simplest example corresponds to a free 
system in a slowly varying chemical potential, E(r, k) = E(k) + V(r), which corresponds to 

H= / dkE(k)a*(k)a(k) + / duV{eu)a*(u)a{u). 



In order to find conditions under which the existence of the limiting / (see ( 2.1 0| ) ) is guaranteed 



we recall the work by [ GMMFj . Suppose the 2-point (which is the only relevant information on the 



state) function of lo z is of the (general) form 

^[a*(x)a(y)} = I d/x(A) V>|(z)VaG/) (3-1) 



where ip\ satisfies the Schrodinger equation (2.11). d[i(\) is an absolutely integrable signed measure. 
On h £ and we assume the same conditions as in ]GMMP | ■ This decomposition includes the closed 



convex hull of pure states (indexed by A). Note, that 

f e (r,k;t) = ( dn(X)w e [iPi](x,k;t) (3.2) 



where w e [ipU(x, k; t) is the Wigner function used by GMMP|] . This relation provides an alternative 



way to derive the Euler equation by quoting their result on w e [-0|](a;, fc; t). 

We write /(r, fc; t) for any limit /(r, /c; t) = lim e |o f e {r, k; i), d s f(r, k; t) = lim e jo d s f £ (r, k; t), s = 
t, r, k and we wish to see what equation is satisfied by such a limiting function. 

Proposition 1 Under the above hypotheses, 

d t f(r, k; t) + 8 k E(r, k)d r f(r, fc; i) = 8 r E(r, k)d k f(r, k; t). (3.3) 



Proof: There are different ways to verify ( |3.3j ). Below we present an explicit and detailed compu- 
tation. 



dtf e (r,k;t) 



-L) / dr ] e^ k d t U- lt [a*{e- 1 r+ ?M £ -V- J)] 



2tt 



-i ( 1 



:»,(/ 



— ( 2^ ) / / d P d P e 



J pipe 1 r-ip'e r+ii (p+p') 
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x J dxdydze- ipx+lp ' v h £ {x 1 z)uol- lt [a*{z)a{y)]~h £ {y 1 z)LO £ e - lt [a*{x)a{z)] 



-(- 

e \2ir 



J drjdpdp' e v 

dydndm e ip ' y - ipm ' lpn uj £ £ - lt [a*(m)a(y)}h £ (n + to, to) 
• J dxdndme- lpx+lp ' m+tp ' n tu £ e - lt [a*(x)a(m)}h E (n + m,m] 



ir)(k-\-u)-\-2ive 1 r 



J dydndm e t{u - v)y - l(u+v){n+m) uj £ ^ t [a*(m)a(y)}h £ (n + m, to) 
-J dxdndme^ u+ ^ x+ ^ u - v ^ n+m ^uj £ £ - lt [a*(x)a(m)}h £ (n + m,m) 

^ / ^ x 2d ,. 

2 d - ( — J / dydndm u) £ -i t [a*(m)a(y)]h E (n + m, to) 



: J dudv 8{k + u)t 



u(n+m-y) 2iv(e- 1 r~ n+ ^ +v ) 



1 

£ I 2^ 

X 



2d 



J dudv S(k + u) 



dxdndm cj £ - H [a* (x)a(m)]h £ (n + to, to) 

(«+m-a) e 2it»(£- 1 r~ "+^+ x ) 



= — y dydndm cu £ - lt [a* (m)a(y)]h £ (n + m,m) 

xe ik(n+m-y) § ( £ -i r _ n + m + y ^ 



+ 



~ (^Y~^J J dxdndmuj £ - lt [a*(x)a(m)]h £ (n + m,m) 

^) / dndydnu e e - H [{a*{y + |)a(y - |)]/i £ (n + y + |,y + 2; 



xe i*(n-Hj)j( e -i r _ " _y) 



~ s\2^) J dnd x dr iu £ £ -i t [{a r {x + -)a(x- -)]h e {n + x- -,x- -) 
xe -ife(«-r/)^( £ -l r _| _ 5 ) 



x/i £ (e r + 



n + n 



n — n 



2 



xh £ (e 



n-n n+n 
1 r+— — ,£ i r — ) 



Now, inserting the definitions of f £ (r,k;t) and h(r,n) into the last expression and expanding 



both at r, we find that 



]imd t f e (r,k;t) 

ej.0 

— ^ lim - [ dndridv e lkr '^ v f e (r -e-,v; t) \e lkn h{r + e 2, -n) - e~ jfc ™ fefr - £-, n) 
27r / eio e J 2 L 2 2 



— ) lim- / d<ndve ikr i- ir > v f e (r,v;t) / f e lfc " /i(r, -n) - e~ ikn h(r, n) 

2n J elO £ J J \ 

+ (2") H ™ y dr]dvi'ne ikr '- ir > v f e (r,v;t)d r J dn e - lkn h(r,n) 



-d k E{r, k)d r f(r, k- 1) + d r E{r 7 k)d k f{r, k; t). 



□ 



4 Remarks 

1. Proposition Jl]) can be extended to time dependent quasi- free time evolutions, U(t, 0) 



-» J" dt' h e (t') 



, where /i £ (i) is a family of self-adjoint operators such that the integral 

gives a well-defined self-adjoint operator. T is the time order ing operation. If h e (x,y;t) = 
h{e(x + y)/2,y — x;et) then the function E{r,k) in equation (3J3) is simply replaced by an 
analogous function E(r, k; t). 



2. The general solution of (3.3) is given by 

f(r,k;t) = f(r(-t),k(-t);0). 

r(t),k(t) are determined by the classical Hamilton's equations of motion: 

d t r = d k E(r, k) 
d t k = -d r E(r,k) 



(4.1) 



(4.2) 
(4.3) 



with Hamilton function E(r, k) and r(0) = r, fc(0) = k. 

3. Analyzing a sum of (properly scaled) quasi-free Hamilton operators is reduced to the sum of 
their respective energy spectra; applicable for example to the case h E — -^{f> + f A(£x)^j + 
eV(ex). 

4. The analysis can easily be extended to generators, H, which are not neccessarily self-adjoint, 
but where the difference between H and its adjoint H* is of order e. As an example consider 
h e = — iA 7 , generating a random walk with bias 7 = (71, . . . , 7^) on the lattice Z d (with unit 
vectors in the positive lattice directions): 



a 

A 7 /(x) = £ [(1 + 7i /2)/(z + <n) + (1 - 1i/2)f(x - (*) - 2f(x) 

i=l 



If 7i = 7(e) be such that 7j(e)/e — > aj as s J, 0. Then, /(r, k; t) satisfies 

dtf(r, k; t) + sin k ■ [(V + a)/(r, fc; t)] = 0. 



(4.4) 



The velocity term comes from the real part of the dispersion relation E{r,k) = 1 — cos k whereas 
the drift a stems from the imaginary part of E. 



6 



5. Notice that the density of particles, p(r, t), is given by 



P (r,t) = / dkf(r,k;t). 



An example where one can say more about p(r, t) is obtained from taking as the initial distri- 
bution the function f(r,k) = w-p{r) with k G [— tt, tt]. Take Z 1 for the lattice case, h = — -|A 
the lattice Laplacian as one particle Hamiltonian and the following initial "diagonal" states cu £ : 

uj E [a*(x)a(y)} = S(x, y) ■ p(ex). 

Then, p(r, t) satisfies the modified Bessel equation 

t 2 d 2 + td t - t 2 A] p(r, t) = 0. (4.5) 

In other words, p(p, t) is the (normal) Bessel function, which for p ^ goes oscillating to as 
\t\ tends to infinity. 



5 Quantum Hard Core Bosons 



In this Section we consider a simple model of interacting Bosons on Z 1 (cf.[ AB]). The pair interaction 
is the so-called on-site hard core interaction (i.e. hard core radius 0), which has the exclusion effect 
that at most one Boson occupies a single lattice site. 

Let -Ha (A) be the following Boson Hamilton operator on an interval AcZ with point interaction 
V{x,y) = A • 8{x,y) of strength A and empty boundary conditions, 

HaW = -\ E b*(x)(Ab)(x) + A £ b*{xfb{xf. (5.1) 

Here, b*(x), b(y) denote the Bose field operators on Z. H\(X) acts on the usual Bosonic Fock space, 
^-"(A). The limit of -Ha (A) as A goes to infinity displays an interesting Hamilton operator. This limit 
is well-defined on states 

Wx)=l[b*(x)\0), Ac A, (5.2) 

xGX 

which define the Fock space, JF /tc (A), of hard core Bosons. |0) is the vacuum state. Let P : J 7 (A) — -> 
!F hc {h.) be the orthogonal projection onto the hard core Boson states. We introduce the following 
annihilation (creation) operators 

a { *\x) =Pb { *\x)P. (5.3) 
In particular, they have the following properties: 

1. The state space ^ rft-c (A) is invariant under b(x) and therefore a(x) = b(x)P. 

2. a*(x),a(y) satisfy mixed (Anti) Commutation relations: 

[a*(x),a(y)}_ = 0, for x ^ y, (5.4) 
[o*(x),o(x)] + = 1, (5.5) 

where [,]-, [,]+ denote the Anti-, commutator, respectively. 

3. lim A ^ 00 i/ A (A) = -±J2 xeA a*(x)(Aa)(x) on T hc {k). 

The commutation relations are familiar from spin systems. Although the previous formulation 
works in any dimension the statistics of the a^(x)'s can only be efficiently "repaired" in one dimen- 
sion. By applying the Klein- Jordan- Wigner transformation (cf. [LSM]) we map the hard core Bose 
field operators onto Fermi field operators. 
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Definition 1 Let 

c(x) = e- l ^<*-i a ' U)aU) a(x), (5.6) 
c*(x) = a*(x)e in ^*-i a * UMj \ (5.7) 

where the sum goes over all j G A left to x. Then, c*(x),c(y) satisfy the usual Anticommutation 
relation. 

Observe that under the above transformation, formally, 

Hx = -\Y J a*(x)(Aa)(x) = -| Yl c*(x)(Ac)(x) (5.8) 
xez xez 

is the free Fermion Hamilton operator. We know that 

n(r, k) = J2 e2lVkc * ( r ~ ^) c ( r + l) ( 5 - 9 ) 

i)£Z 

are the locally conserved quantities for the free force evolution. Rewriting them in terms of hard 
core Bosons we get 

„(r, k) = 53 e 2lvk e N ^- 71 - r+7l) a*(r ~ 7 ? )a(r + rj), (5.10) 

r/eZ 

where N(x, y) = Yl x <j<y a *(i) a (j)- Although the time evolution of the a*(a;)'s is quite complicated 
(and unknown) we are only interested in a special combination appearing in n{r, k). The problem is 
then simply reduced to free Fermions on the lattice (but only as far as this local quantity, n(r, fc), is 
concerned). It is then a simple matter to apply Proposition 1 : 

Proposition 2 Suppose that f(r,k) = lim e | w e [n(e _1 r, k)] exists (for example if the functions 
in the decomposition of uj e (see \3.\ ) are e-oscillatory (see fflMMl\ ) and bounded). Then, under the 
dynamics generated by the Hamilton operator \5.f\ ), f(r, k; t) exists and satisfies the free force Euler 
equation 

d t f(r,k;t) + sinfc • d r f{r,k;t) = (5.11) 

Remarks: 

1. Since we are working in an infinite discrete configuration space we can map any configuration 
of Bosonic hard rods with positive (integral) length bijectively onto a configuration of Bosons 
with core length 0, i.e. the above case. Therefore the last proposition extends to hard core 
interactions of positive length. 

2. In continuous configuration space things are much more delicate. Whereas the same analysis 
is literally applicable to Bosons on R 1 with core radius 0, the hard rod model with Neumann 
boundary conditions (for example) seems to be very different and rather complicated. 



6 Fluctuations of the Density 

The Euler equation corresponds mathematically to a law of large numbers. The simplest equation 



for the density profile is (4.1:). Now we study (non-equilibrium) fluctuations of the particle density 
for the quasi-free systems of Section 3. The field of fluctuations, £, £ (ip, A, t), of a local observable, A, 
is defined by 

FU>, A, t) = e d ' 2 J dr i>(er) [r r ^- H (A) ~ ^(t^- h (A)] . (6.1) 

Here, tp is any test function and uj e a family of states as in Section 2. Further, we denote by T r j(A) 
the space and time shift of the observable, A, by r and t, respectively. We are interested in the 
density, r r ,t(A) = at(a*(r)a(r)), in particular in realizing the limit 
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hmf {ip,A,t) = I drij;(r)£(r,t). 



(6.2) 



We should therefore study the characteristic function, uj e (e 1 ^ vl'^vY and try to reconstruct the 
limiting dynamics of the fluctuation field. Since we restrict our attention to quasi-free states 
(cf.[BRj) and we always assume nice clustering properties it is fair to expect Gaussian behavior. 
The precise meaning of this in the quantum case (the so called quantum central limit theorem) 
can be found in |GVVl| and |GVV2fl . We are satisfied therefore with studying the covariance 
lim £ 4o w e [£ e (V>, t)£ e ((f), s)] with £ £ ((f),s) given by (6.1) for A = a*(0)a(0) and for a class of test func- 
tions ip,(f> £T>. 

Clearly, the dynamics of the fluctuations can only be derived when the law of large number is 
already established. Therefore, while the next proposition does not follow directly from Proposition 
1, it is essential to make sense of fluctuations that one first understands the averaged behavior. 
We put ourselves therefore in the same context as for Proposition 1 and we have that f(r,k;t) = 
lim e |o f £ (r, k; t). On the level of the one-particle dynamics we recall that we denote by r e (t) the 
rescaled position operator, see (2.3), and we assume that it converges to an operator r(t) as e \ 0, 
SCO pMP| , @, |S§. This convergence takes place on a suitable domain V dense in L 2 (R d , dx). Let 
us write 

I \d/2 



(e ir ^)(x,k) = (j^j Jdqi>(q)< 



iqr(x ,k;t) 



(6.3) 



where r(x,k;t) is the solution of (4.2). (If the motion is force free then (e lr (')?/>) (x, k) = ip{x + 
d k E{k)t).) 

Proposition 3 Under the above conditions on the initial states and on the Hamilton operator, 

limw £ [f (il>,t)£ e (<M] = / dxdk (j r(t ~ s) i>)(x,k)ct>{x)f(x,k;s)(l±f(x,k-,s)), (6.4) 
where +/— concerns Fermi/Bose statistics, respectively. 



Proof: 



u> e rw>,*)r(<M] 



drdr' 



' ip(er)<j)(er'){uj £ [a*(r, e~H)a(r, e -1 i)a*(r' ,e' 1 s)a(r' , s^s)} 
N (r, e-H)a(r, e" 1 *)] u; £ [a*(r', £- 1 s)a(r', e _1 s)]} 
e d I drdr'i/;(er)0(£/)w e [a*(r,£- 1 i)a(r',e- 1 s)]w £ [a(r,e- 1 i)a*(r',e- 1 s)]. 



By changing the order of a and a* in the last expectation we get two parts. The one including 
the (anti) commutator is equal to 



e d I drdr' i)(Er)<t>(er')u s [a*(r,e-H)a(r' ^^s)]^ ([a(r, £~H), a*(r', e -1 *)]±) 

^_L^ 1 £ d J dr'dq^{q)(t>{Er') J dzdr e lqer e -i^(t-«)« _I ( r , z ) ^M*-^" 1 ^ 
XLu £ [a*(z, £~ 1 s)a(r', £ _1 s)] 



2tt 



d/2 



d/2 



dr'dq^{q)cj>{er') / dz e 1 ^^'^ {r' , z) uo^^a* {z)a{r')\ 



i-] ' e d I dxdr,dq$(qMe(z - \))e***«-'\x -\,x+ f K_ ls [a*(z + ^)a(x |) 



2tt 



d/2 



dxdr]dkdqij(q)(t)(x - ehe**"^ (e~ l x - ^e^x + |) e - ir > k f e {x, k; s) 



9 



Now remember our scaling of the matrix elements of h e and take the limit e J. of the last 
expression: 



i-) J dxdkdq i>{q)4>(x) (J ' dr, e iqr ^ (x, 7?) e"**) /(z, fc; a) 
(e^'-ty) (x,k)cj>(x)f(x,k;s). 



dxdk 



For the second part (forgetting the (27r) d ^ 2 ) we need 



lime d y dqdr' ijj(q)<j)(er') J dr e iqer u e s -i s [a*(r, e' 1 ^ - s))a(r')] u e s - ls [a*(r')a(r, e- l {t - a))] 

= Kme d [ dqdr' i>(q)(t>(er') J drdzdz' e lq£r e-' he ' 1{t - s) (r,z)u; £ E - ls [a*(z)a(r')} 

xe^O" 1 V ",r)u E £ - ls [a*(r')a(z')} 
= lim^y dqdr' i>(q)ct>{sr') J dzdz'e i « £r( ^ 1(t - ;i)) (z , ,z)^_ ls [a*(z)a(r')]<- 1;5 [a*(r')a(z')] 
= hm £ d / dg<fc*/#?)(«e(x - |)) | dz> e^^Hz' , x + |) u^ ls [a*(x + |)o(x - |)] 

x W |_ ls [a*(x - |)a(*')] 
= lim£ d / dqdxdri4>(q)(t>(ex) [ dz'dk e lqrS{t - s) (z' , x + -) e~ lvk f e (ex, k; s) 

xwf_ ls [a*(a;- -)a(/)] 

= lime d / dqdydr]df]dkdk^(q)(j)(e(y + h) e iqrE {t ~ s) (y -%y + 
e|0 7 2 2 2 

- Um / dqdydr]df]dkdk^(q)<f>{y + £ h e^^ie^iy + e^-) - "^fi , e- 1 [y + e^-) + V + ^ * 



eioj * " ' ' 2 y 2 / 2 ' v " 2' 2 ' 

xe-^-^rd/ + e^), k; S )f(y, k; s) 

= Um / dqdydr]df]dkdk^{q)(p{y + e^-) e lqr{t - s) (y + e^-, 1 ^^-) 
ej.0 7 2 2 2 

xe-*»*-^*/ e (i/ + fc; *)/ £ (y, *; *) 

= lim f dg^di/dPdpdp^^^ 
eJ.0 7 

= (27r) d/2 y d^dpdp (e^'-'ty) (2/,p)0(y)e- lpp 7(y,p + p;s)/(y,P~p; S ) 
7 Particles in a Periodic Potential 

So far we have not dealt with the important case of particles moving in a periodic potential. More 
precisely, where the system is quasi- free with one particle Hamiltonian of the form h = — A + U on 



R d with periodic potential U. Recently, Spohn (|S3|) studied their long time behavior. To obtain 



the corresponding Euler equations in this case is clearly of interest but some problems immediately 
arise. 
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The first one is related to the identification and the proper interpretation of the corresponding 



Wigner distribution function. In other words, by what to replace the equations (2.7)-(2.9)? We do 
not know a unique good answer to that question (but see the discussion in pMMFfl , TMMP| , |§ 



and AKQ. An obvious generalization of (2.4) goes as follows. 

Let ipnk be the eigenfunctions of the Hamiltonian H = — A + U, with periodic potential U, having 
eigenvalues E n (k). Introduce the generalized Fourier transformation (cf.flRSf) 



X(k,n) = / dxi(j nk (x)x(x). (7.1) 

We consider now an ensemble (we denote creation and annihilation operators by a*(x),a(y)) of 
particles subject to a periodic potential U and a family of states u £ for which the following limit 
exists (A* is the first Brillouin zone) 

f n (r,k) =lime / dv e- irv uo e [a*{k + ej, n)a(k - e~, n)], rGR^eA* (7.2) 

elO J R d 2 2 



and call it the n-th band distribution function of the ensemble. d(k,n) is to be understood as (7.1), 
but now as operators. We recover (2.7) at time t = if in (7.1)-(7.2) we use tp n k(x) = 



Comparing with (2.8)-(2.9), we have 



7k. r 



and 



dkf n (r,k)=limui e (a*(e 1 r,n)a{e 1 r 1 n)) (7-3) 

e|0 



drf n (r,k) = lira eu E (a* (k , n)a(k, n)) (7-4) 

elO 



Proposition 4 Let at(-) be the time evolution of the ensemble generated by the one particle Hamilton 
operator H = —A + U on R d : Hipnk = E n (k)ip n k ■ Then, 



f n (r,k;t) =lime / dv e~ irv {uj e o a e -i t )[o*(fc + e|, n)a(k - e J, n)} reR^eA* (7.5) 



elO J Rd v " WL v 2' ' x 2' 

exists and 

d t fn(r, k; t) + d k E n (k) ■ d r f n (r, k; t) = 0. (7.6) 
Proof: This is an easy task since ata(k, n) = e~ ltEn ^d(k, n) and/ n (r,fc;t) = f n (r — dkE n (k)t,k).D 
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